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Abstract 

The gravitational neutrino oscillation problem is studied by con- 
sidering the Dirac Hamiltonian in a Riemann-Cartan space-time and 
calculating the dynamical phase. Torsion contributions which depend 
on the spin direction of the mass eigenstates are found. These effects 
are of the order of Planck scales. 
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1 Introduction 



Neutrinos are the only particles for which a poem has been written; " Cosmic 
Gall" by John Updike, quoted in Recently neutrinos have attracted a 
lot of attention in high energy physics; see ^ ^ and references therein. A 
major problem at present is the problem of solar neutrinos. Thermonuclear 
reactions in the sun produce large number of electron neutrinos; the standard 
solar model allows us to predict the flux of these neutrinos approximately 
as 7.3 ±2.3 SNUQ. On the other hand experiments which have been done 
for decades give a measured flux of 2.55 ± 0.17 ± 0.18 SNU - approximately 
one third of the theoretical rate. This is the solar neutrino problem. A 
much-publicised solution is "neutrino oscillation (or mixing)". According to 
neutrino oscillation, in contrast to the conventional view that neutrinos are 
massless particles, they actually have a mass different from zero, although 
very small. The essential idea belongs to B.Pontecorvo 0. Technically 
speaking it is that the neutrino states of definite mass do not coincide with the 
weak interaction eigenstates. In other words, an electron neutrino emitted by 
sun becomes a linear combination of all three neutrinos (electron, muon and 
tau) while travelling in space towards the earth. Therefore the probability 
of detecting z/g emitted by sun as t'^ or Ur on the earth is different from 
zero. This situation, of course, invokes physics beyond the standard model. 
Although there are three neutrino families, the essential physics of neutrino 
oscillations is illustrated very well by considering the interactions of only two 
of these, which we choose to be I'e and f^. 

All the above arguments have been cast in Minkowski spacetime. However 
^SNU stands for solar neutrino unit and 1 SNU = 1 capturc/scc/lO"^^ target atoms 
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we know that we live in a curved spacetime - perhaps even in a curved 
spacetime with torsion. Therefore, in more recent years, physicists have 
turned their attention to specifically gravitational contributions to neutrino 
oscillations - see |^ H, ^ ^ ^ and references therein. 

In this paper we reinvestigate the effects of gravitation, but we consider 
a spacetime with torsion, as well as simply curvature; this is the Einstein- 
Cartan theory, or more specifically the Einstein-Cartan-Dirac theory in our 
case, since we will be considering Dirac particles in the spacetime. Some 
work along these lines has already been done |ll|], but our results do 
not entirely replicate these findings. The essence of our work is to calculate 
the dynamical phase of neutrinos, by finding the form of the Hamiltonian, 
H, from the Dirac equation in Riemann-Cartan spacetime. The phase then 
follows from the formula 



In what follows, -0 is a Dirac spinor and is a 4 x 4 matrix; the expo- 
nential function is therefore well-defined, as usual, by the series expansion. 
The Hamiltonian H will depend, for example, on momentum p, and this is 
expressed not as a differential operator but simply as a vector. 

After we have introduced briefly the Einstein-Cartan-Dirac theory in 
Sec. 2 we flnd the hamiltonian of a Dirac particle in subsec.2.1. In sec. 3 
we investigate, in turn, the neutrino oscillations for azimuthal and radial 
motions. The conclusions are given in sec. 4. 




^(0) . 



(1) 
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2 Einstein-Cartan-Dirac Theory 

A Riemann-Cartan space-time is defined by the triple {M,g, V} where M is 
a 4- dimensional differentiable manifold, equipped with a Lorentzian metric 
g and a metric compatible connection V. The metric tensor may be written 
in terms of orthonormal basis 1-forms {e"} as 

g = Vabe"" = -e° e° + ® + ® + ® , (2) 

while the connection is specified by a set of connection 1-forms {T^f,}- Metric 
compatibility requires 

Dr]ab = dr]ab - ^"aVcb - ^"bVac = , (3) 

so that Tab = —Tba- The Cartan structure equations below define the torsion 
2-forms {T"} and the curvature 2-forms {R^b} of the space-time: 

De" = de" + T\ A = T\ (4) 
DT\ = dT\ + T\AT\ = R\. (5) 

In this study d, D, i^, * denote the exterior derivative, the covariant exterior 
derivative, the interior derivative and the Hodge star operator, respectively. 
The local orthonormal frame {Xa} is dual to the coframe {e"}; 

tx^e' ^ z.e" = e\Xa) = 5^ . (6) 

The space-time orientation is set by the choice eoi23 = +1- In addition, the 



Dirac matrices 7" are in the standard Bjorken and Drell form |jT2| but for 
our metric convention they satisfy 

7^7' + = -2r7"''l4x4 . (7) 



We use the formalism of Clifford algebra-valued exterior forms for the 
Dirac particles. In this formalism 7 = 7^60 and *7 = 7*^ * Cq. Furthermore, 
the covariant exterior derivative of spinor fields in the Riemann-Cartan space- 
time is 

The connection 1-forms, {F";,}, are decomposed as follows 

r\ = u\ + K\ . (9) 
Here, {uj°'b} are Levi-Civita connection 1-forms which are given by 

de" + A e'' = (10) 
and {K"'i,} are kontorsion 1-forms which are given by 

K\ A = T'^ . (11) 
The Einstein-Cartan-Dirac Lagrangian density 4-form is 
L = -^Rab A *(e" A e^) + ^{^ * 7 A + A *7V'} + ^^V' • (12) 
The field equations are then |Q , 0] 

ih* -J A {Dtp -^Ti/j) + *mci/j = 0, (13) 



I A *{ea A Cfe A = — A Cf, A ^^775^/^ , (15) 

where T = laT"" is the torsion trace 1-form and £ = \J~^ ~ 10"'^^ m is 
the Planck length. As is seen from eqn. ([T5|) the source of torsion is spin and 



magnitude of torsion is proportional to square of the Planck length. In this 
work, the spinorial contributions to the curvature will be assumed small, and 
therefore neglected. That is both the Dirac stress-energy tensor on the right 



hand side of eqn.(|T^) and the torsion dependent terms on the left hand side 
are dropped. 

The Dirac equation (^) can be recast as 

ih-f^'iaiD^ - ^Ttp) - mctp = . (16) 

Here the torsion trace 1-form T is cancelled by a similar term in the covariant 
exterior derivative and what is left behind are the totally anti-symmetric axial 
components of the torsion. These are going to be treated as given background 
fields. 

2.1 Hamiltonian of a Dirac Particle 

As mentioned above, our aim is to calculate the phase of a neutrino beam 
in Riemann-Cartan spacetime. We do this by finding the Hamiltonian for 
a Dirac particle. The spacetime we consider is Schwarzschild spacetime and 
the metric is therefore given by 

g = -fdt ®dt + ^dr ® dr + r'^^dO ® d9 + sin^ Od^p ® di^) . (17) 
/ 

The orthonormal basis coframe is (from (H)) 

= fcdt = jdr = rdd e^ = rsm6d(f (18) 



where = 1 — Using (|TB|) and (jT^) we can calculate the Levi-Civita 



connection 1-forms, uj°-f,. We then substitute eqn.(^ and the Levi-Civita 



1-forms into eqn.(^) and then into eqn.(|T^). Now we use the identity 

7I7', 7l = + + 22e'^'"=S57d (19) 

and eqn.(|TTl), and then write the contorsion 1-forms in component form as 
K"'b = K°'bcG'^- The resulting Dirac equation is 

Ah n d .^r 1 f, ih cot 6. 

''^ Y'l^-mc+7AW'-'7574i/' = (20) 



r sin 6 dip 4 

where 75 = —ij'^'y^'y'^'y^. From the last term we see that it is only the totally 
anti-symmetric part of the contorsion tensor (in other words, the purely axial 
component) which contributes to the Dirac equation. [Actually, in theories 
in which spin is the source of torsion, the contorsion tensor is in any case 
totally anti-symmetric ll^.] Relative to the orthonormal basis 1-forms {e"}, 



we may write 

Kabc = eabcfA^ (21) 

for some vector fieldQ A = A°-Xa- Noting that eabcf^"''^"^ = ~3!5j and defining 
the components of the momentum operator by 

d 1 f 

Pr = -^% + - + ^) (22) 
ih , d cot 9, 

r sm U 0(p 



2 From cqn.(|l8|) Aq = j^At , Ai = fAr , A2 = ^Ag , A3 = t^^.^- This means 
A — {Ai, A2, A3) are not the cartesian components, rather the spherical components in 
the orthonormal basis. 
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the eqn.(^) becomes 



+fmc' + ^hcf^.^aA''}^ . (25) 





d = 


(f 


:) 








and S = . Comparing the above with the Schrodinger equation 



= H^l^ (26) 



we find that 



H = fa'prc + fa^pec + fa^p^c + f fimc^ + ^^c/A°75 + ^^c/S.A . (27) 



3 Neutrino Oscillations 

Neutrinos that are produced in the sun arrive on the earth. According to the 
theory of neutrino oscillations, a neutrino with fixed fiavour which is emitted 
by sun becomes a superposition (in our model) of two neutrinos, i.e., they 
mix during the journey. In this case, even if only one kind of neutrino were 
emitted by the sun, the probability of measuring the other kind of neutrino 
on the earth would not be zero. This is called "neutrino oscillation" and is 
measured by the interference pattern observed on the earth. The starting 
point of neutrino oscillation theory is that fiavour eigenstates are different 
from mass eigenstates; each flavour eigenstate is a linear superposition of 
mass eigenstates; 

Ve = cos 6z/i + sin 6z/2 (28) 
= — sin Gz/i + cos Qi'2 (29) 



where Pe and denotes electron and muon states, respectively, and fi and 
V2 mass eigenstates, and 6 is a mixing angle. It should be clear that this 
scheme only works if the neutrino masses are different from each other, and 
therefore in general are non-zero; this means that in general there are right- 
handed neutrinos as well as left-handed ones. The right-handed neutrinos, 
however, interact with matter only very weakly. That is why they have not 
yet been observed. 

Now how does the interference work? The wave packets which represent 
electron and muon neutrinos are constituted by the sum of the wave packets 
which represent mass eigenstates. These wave packets will travel with dif- 
ferent speeds, hence, even though they may be centered at the same point 
in space at the point of production, there will be a seperation between the 
packets at the point of detection. In order for interference to be observed, 
however, the spatial separation between them must not be too large. The 
waves will then interfere at the detector and we can observe the phase dif- 
ference. For more detailed discussion see Refs. 0, |T^ . 

The model which explains this phenomenon realistically is the radial mo- 
tion analysis. We begin the investigation, however, by studying azimuthal 
motion. Although not realistic, this has the virtue of simplicity; and it has 
also been considered in Ref . [|l^ , though in a different formalism. Our method 
of procedure will be to write down the Dirac equation and find phases cor- 
responding to mass eigenstates, then finally calculate phase differences. 
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3.1 The Azimuthal Motion 

The Hamiltonian for the azimuthal motion, that is, p= {pr,Pe,Pifi) = (0, 0,p) 
could be written from eqn.(p7D as 

H = fa\c + f(3mc^ + ^/ic/A%5 + ^^c/S.A . (30) 

This Hamiltonian couples the positive energy states to the negative energy 
states because of and 75 : 

H'4) = Eijj (31) 

where is a four- component spinor whose first two components correspond 
to positive energy and second two to negative energy states. Firstly we 
decouple these states (or block diagonalize the Hamiltonian) by writing 

det(H-E)=0 (32) 

E+ 

where E = \ I . Therefore the upper block E^ and lower block 



E^ 



E^ satisfy 



E^ = fpc+l^ + hcf{-A'a' + a.A) (33) 
E_ = -fpc-L_^ + ^ncf{-A'a' + a.A). (34) 

Here we assume that the torsional terms are small compared with the mo- 
mentum and mass termsQ. With a unitary transformation 

ij = U^ (35) 



Numerical values are given at the end of next subsection 
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where U is a unitary matrix and ^ = | | , our equation system decouples 
into two equations 

= (36) 
E-^- = th^^^_. (37) 

We will be interested in only the positive energy states for simplicity. Similar 
calculations may be done for the negative energy states. At the end these 
states will correspond to the mass eigenstates whose linear superpositions will 



define, for example, electron and muon neutrinos as given by (p8|) and (pOl). 
Now eqn.(|36D can be cast easily as two separate equations by diagonalizing, 
one for spin up, one for spin downQ; 

{fpc+^ + hcfAUl = ^hlcl (38) 

{fpc+^-lhcfA}^^ = (39) 

whose time evolutions are given by 

d(t) = e-*^W^l(0) (40) 
Ciit) = e-'^'mliO) (41) 



where A = y {A^ — A^y + {A^)^ + (^4^)^. The phase of the spin up state is 
given by 

^\t) = -J {fpc + L-— + %cfA}dt . (42) 



^ Spin operator is chosen as S. Besides, it is only in the MASSLESS limit that the 
neutrino helicity is ±1, so the state we are discussing is actually allowed, since we are 
discussing MASSIVE Dirac neutrinos. 
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For ultrarelativistic neutrinos pc c::^ E and cdt ~ Rdip , so the above becomes 
IJER^ fm'^c^R^ 3,,,^, , 

Similarly the phase of the spin down state is 

$^ = ^{^A^ + ^^^Av. - IhfARA^} . (44) 

These phases alone do not have an absolute meaning; the quantities relevant 
for the interference pattern at the observation point are discussed below. 

Now consider electron and muon neutrinos given by ( pSf ) and ( p9D above, 
where ui and z/2 correspond to ^+'s. If electron neutrinos are produced at 
t = 0, the time evolution of Ue is 

fe(^) = cos 9z^i(t) + sin Oz/2(^) 

= cos Oe-'^'^'Ki + sin Ge^^^^^^Va . (45) 



The probability of measuring muon neutrino at a later time is given by 

Pii^e T^ii) = sin^ 29 sin^ ^ (46) 

where A$ = $2 — As is seen from (|^) the phase difference has an 
absolute meaning. Now there are four possibilities. 

• Both mass eigenstates are spin up, in which case eqn.(^5D becomes 

i/e (t ) = cos ee-**^ 1^1 + sin Ge'^^^ . (47) 
From (|43| ) and (^If), the phase difference is 

A$ = $^ - $1 = ^^i?Av. . (48) 
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Both mass eigenstates are spin down. In this case an analogous equa- 
tion to (^7|) holds and the phase difference A$ is 

A$ = $i - $1 = ^^i?Av. . (49) 

The first mass eigenstate is spin up, the second spin down, then 

A$ = $^ - $1 = - 3fA}RA^ . (50) 

The first one is spin down, the second spin up, then 
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A$ = $^ - $1 = + 3fA}RAcp (51) 



where Am^ = ml — mf. From (58) and ( ^91) it is seen that if both mass 
eigenstates have the same spin projection there is no contribution to the 
oscillation coming from the torsion, but if the mass eigenstates have opposite 
spin it is seen from (^0]) and ( |51| ) that there is a contribution to the neutrino 
oscillation coming from the torsion. 

When we do an order of magnitude analysis we see that the contribution 
coming from the torsion is really very small compared with that coming from 
the mass difference^. Here we also emphasise that j is the locally measured 
energy 0. If there is no torsion, i.e., A°- = 0, our result becomes the same as 
the result of Ref . , and in this case there is no gravitational contribution 
to the phase difference, since the locally measured energy is j rather than 
E. We show below that in the case of radial motion the situation is different. 



Numerical values are given at the end of the section 3.2 
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3.2 The Radial Motion 

The Hamiltonian for the radial motion, that is, p = (p, 0,0) may be written 
from eqn.(p7D as follows 

H = fa'prc + fPmc' + hcfA^j, + hcf^.A . (52) 

Applying the same arguments as in the previous subsection we find the upper 
and lower blocks of the Hamiltonian 

= fpc + ^ + ^hcf{-A'a' + a.A) (53) 
2p 2 

E- = -fpc-^-^hcfi-A'a' + a.A). (54) 

As before, we are only interested in the positive energy states which obey 

= ih^^^+. (55) 

Again we first diagonalize eqn.(^) in order to find eigenvalues of spin up 
and down states. In this case the phases appropriate to the spin up and spin 
down particles are, respectively, 

. IK, 2MGE ^ tb m^c^ ^ 
$T = -|— Ar — In \ Ar 

+ -M(Ar-^ln^)} (56) 
2 ta 

$i = -{— Ar — In— H Ar 

he & ta '^E 

--M(Ar-^ln^)} 57 
2 va 



where A = \J (A° — A^)^ + (^^)^ + (^'^)^- Here we made the assumptions 
pc ~ E, cdt ~ dr, / = (1 — ^^)^/^ ^ 1 — Ar = Tb — ta, where A and 
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B refer respectively to the points of production and detection; hence is 
the radius of the sun and is the distance from the center of the sun to the 
surface of the earth. 

If we apply the same procedure as in the previous subsection, from 
eqn. (|56|) and (|57D there are three different phase differences; 

• both mass eigenstates have the same spin state 

A$ = $^ - $1 = $^ - $1 = ^^Ar (58) 

• the first is spin up, the second spin down 

A$ = $^ - $1 = ^!^Ar - 3A(Ar - ^ In ^) (59) 



• the first is down, the second up 

A$ = $^ - ${ = Ar + 3A(Ar - ^ In ^) (60) 

where the mass difference is Am^ = m\ — ml. 

If an order of magnitude analysis is done, for vacuum Am^c^ ~ 10"^" eV'^, 
E 10 MeV, ta = Rq = 7 X 10^ m, the earth-sun distance Ar ~ 
1.5 xlO"*^^ m and = Ar+r^ ~ 1.5x 10"*^"^ m and for the sun ~ 1.5 km. 
As for the magnitude of the vector A^, it follows from (|15D that 

A"" ~ fii^ij) . (61) 

Then ~ 10^^^ m~^ and hence the torsional term in (^), ( |60D is much 

smaller than the mass difference term. When we put ^4" = we recover the 



same result as Ref.|TO|, but we see that in general there is a contribution to 



neutrino oscillation coming from the torsion of spacetime. 
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4 Result 



We have treated the gravitational neutrino oscillation problem in a different 
way from previous authors by calculating the Hamiltonian of a Dirac particle 
in Riemann-Cartan spacetime and finding the dynamical phase. We began 
for simplicity by treating the (unrealistic) case of azimuthal motion; this 
gave us the opportunity to compare with in which azimuthal motion 
is considered in Riemann spacetime. We have found that there is a tor- 
sional contribution which depends on the spin directions of mass eigenstates, 
eqn.(|50|) and (|5lD. If we put the torsion equal to zero in these equations we 
recover the same result as Ref. ||1U[|. Next we studied the radial motion which 



is a more realistic model for our problem. We found here also a torsional 
contribution to the neutrino oscillation that is dependent on the spin polar- 
izations of the mass eigenstates. Again, when the torsion vanishes we recover 



exactly the same result as in ref||TO| 
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